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We study surface states of ID photonic crystals using a semiclassical coupled wave theory. Both 
TE and TM modes are treated. We derive analytic approximations that clarify the systematics of 
the dispersion relations, and the roles of the various parameters defining the crystal. 
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I. INTRODUCTION 

Photonic crystals are artificial low-loss dielectric struc- 
tures with periodic modulation of refractive index, which 
have attracted considerable attention in the last two 
decades. Due to Bragg reflection, electromagnetic (op- 
tical) waves cannot propagate through such structures 
in certain directions, at certain frequencies. Hence, pho- 
tonic crystals can control the propagation of electromag- 
netic waves in novel ways, with obvious application to 
dielectric mirrors, dielectric waveguides, and dielectric 
laser cavities. 

As a way to efficiently inject light into a photonic crys- 
tal (PC) waveguide, it has recently been proposed to use 
surface electromagnetic waves (SEW) 00. In those pa- 
pers, the photonic crystal was a two dimensional array 
of rods, of infinite length normal to the plane of inci- 
dence. Instead, we have studied SEW on a semi-infinite 
one-dimensional (ID) photonic crystal sketched in Fig. 
1. While retaining all useful properties of 2D and 3D 
photonic crystals, a ID dielectric periodic structure with 
high refractive index contrast is more attractive from a 
technological point of view. 

The usual theoretical methods for wave propagation in 
ID photonic crystals, including SEW , are the Floquet- 
Bloch modal formalism, coupled wave theory, and the 
transfer matrix method. Among these three, the coupled 
wave approach^, 0,0 offers superior physical insight and 
gives simple analytical results in limiting cases. Unfortu- 
nately, the conventional coupled wave theory of Kogclnik 
fails in the case of high refractive index contrast, which 
is essential for a functional ID photonic crystal. 

In this paper, we apply our recently developed semi- 
classical version of coupled wave theory^, 0, to SEW 
on ID photonic crystals. The method is analytically al- 
most as simple as conventional coupled wave theory, and 
is practically exact for the achievable ratios (e.g. 1.5:4.6) 



of the indices of refraction of the materials available to 
build devices. We present here a unified description of 
TE and TM SEW. A detailed account of the properties 
of the TE surface modes has recently been given by us 
in Ref . jfj] ; here we complement these findings with those 
for TM modes, which are slightly more complex due to 
the presence of Brewster points in the bandgaps. As a 
result, we thoroughly clarify the systematics of solutions 
for surface EM waves in semi-infinite ID photonic crys- 
tals. 

Our method is formally quite different from that re- 
cently presented in Ref. or those in Ref. [ll|> so 
in Section II we provide a short summary of the transfer 
matrix approach, in the notation of our previous work0. 
In Section III we rederive the exact equations for SEW 
of TM modes and obtain from them various analytic ap- 
proximations for a semi-infinite crystal. The analogous 
equations for TE modes were given in Ref. . With these 
in hand, we discuss systematics of SEW. In Section IV we 
apply the semiclassical approximations of Refs. and 
to surface waves, and show that the second approxima- 
tion is very accurate both for the dispersion relation and 
the bandgap boundaries. 



II. TRANSFER MATRIX METHOD FOR A 
PERIODIC CRYSTAL 



We wish to describe surface states that form at the 
interface between a medium of low refractive index, no, 
and a semi-infinite 1-D photonic crystal with layers of 
refractive indices ri\ and ri2 and thicknesses d\ and di. 
We choose a coordinate system in which the layers have 
normal vector along OZ. As shown in Fig. the crystal 
is capped by a layer of the same material but different 
width, d c . For convenience of presentation, we split this 
termination layer of index of refraction n\ and width d c 
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into two sublayers, of lengths d s + d t = d c . The first 
sublayer extends from z = — d s to z = 0. Then the 
periodic array that forms the ID photonic crystal consists 
of "cells" each made of three uniform layers of widths d t , 
g?2 and d\ — dt whose respective indices of refraction are 
n\ , ri2 and n\. (If dt = d\/2, the unit cell will have 
reflection symmetry, which simplifies some relations, but 
does not change any physical results.) The first cell, given 
index 0, ranges from z — to z = d = d± + d2', the second 
is given index 1, and ranges from z = d to 2d, etc. The 
p-th cell runs from pd to (p+ l)d and has n(z) — rv\ when 
pd < z < pd + dt or pd + d t + < z < (p + l)d and 
n(z) = ni when pd + d t < z < pd + d t + d2- We choose 
m > ri2 > no. 

For monochromatic TE waves the electric field is par- 
allel to the OY axis. As in Ref. 0, we write 



E = E y (z)e i 



J(kj3x — tot) 



(1) 



and 
A™ 



iki di 



COS fc 2 G?2 



n\k 2 



2 \n\k\ n\k 

gTM _ ^_ ik 1 {d 1 -2d t ) I "\ h ~ 



nkk\ 



sin k 2 d2 



njk 2 

n\k\ n\k2 



sin fc2G?2 



(7) 



Once these are known, the Bloch waves of the ID crystal 
are determined by the eigenvalue equation 



[ a p+i 



A 1 z 



M 1 z 



and therefore the A satisfy 

(A- X)(A* - X)-\B\ 2 = 0. 
Using that det M = \A\ 2 - \B\ 2 = 1, one finds 



(8) 



(9) 



where to is the angular frequency, k — lo/c. is the vac- 
uum wavenumber and (5k is the (constant) ^-component 
of the wavevector of modulus k(z) = Lun(z)/c. For an 
electromagnetic wave entering the ID photonic crystal 
from a uniform medium, one has 



= no sin 8q, 



(2) 



where 9q is the angle of incidence measured from the 
normal. For monochromatic TM waves it is the magnetic 
field which lies parallel to the OY axis. Following Ref. 
we write 



H = H y (z)e y e 



i{k(3x — tot) 



(3) 



For piecewise constant e(z) = n 2 {z), the solutions of 
Maxwell's equations for E(z) and H(z) in the p-th cell 



are 



Ey(z),H y (z) = a p e lkl{z - pd) +b p e- lkl{z - pd) 



ik\ (z—pd) 



+ f P e 



-ik± (z—pd) 



(4) 



in regions pd < z < pd + dt, pd + dt < z < pd + dt + d 2 
and pd + d t + d 2 < z < (p + l)d respectively. Of course, 
the coefficients (a p , b p ), (c p , d p ), and (e p , f p ) differ for 
TE and TM waves. 

The transfer matrix, M, is defined so that 



a p +i \ _ 



b p+ i 



= M 



p = 



A B 

B* A* 



(5) 



Matching the fields in eq. (@J and their derivatives (the 
derivatives divided by n 2 ) for TE (TM) waves, one finds 



A± = Re(A) ± v/Re(A) 2 - 1, 
and eigenvectors 



B 

X+-A 



(10) 



(11) 



In the bandgaps A± are real since \ReA\ > 1. In contrast, 



in allowed bands A± 



±id 



with 4> real. The bandgap 



boundaries are at ReA = ±1. Furthermore, from eq. 
(I10JI A_A + = 1. For surface states one chooses the A that 
corresponds to a damped Bloch wave when z — > +oo. It 
must fulfill the condition |A| < 1, and is A_ ( A + ) when 
Re(.A) > 1 ( < —1.) For simplicity we write it simply as 
A from here on. 
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FIG. 1: Sketch of a typical ID PC, for the case used in 
calculations (ni = 2.0, di — 100 nm, ri2 = 1.5, c?2 = 250 nm,) 
A symmetric unit cell was chosen, and d a = 80 nm. 
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Brewster Points 



Brewster points exist only for TM waves, and require 
Re A™ = ±1. From eq. we find that 



Re A = cos kidi cos k 2 d 2 
1 / n\k 2 n^k\ 



2 \n 2 k\ n\k 2 



sin£;ie?i sin fc24- (12) 



It can be easily checked that a Brewster point occurs 
when 

= — i and k\d\ + k 2 d 2 = ran, (13) 

with m = 1, 2, ... an integer which we assume to be the 
bandgap index. The first of these equations determines 
[3. In particular, we have 



y/n{ - P 2 _ V^hT 2 



and, as a result, 



nin 2 



(14) 



(15) 



The second of equations (|13|1 then fixes the value of k to 
be 



k Br = 



d ly /nl-^ + d 2 ^n 2 2 -p 



(16) 



For our reference case, see Fig. ^ we nn d that Pbt = 
1.2, ksr — 0.00816 nmT 1 when m — 1 (first bandgap). 
A numerical determination of the bandgap boundaries 
confirms that the gap width shrinks to zero at this point. 



note that the ratio qo/k% depends only on /3, and is in- 
dependent of k. Furthermore, from eq. 10 one has that 
A — A, B and B, have the same modulus and therefore 
eq. (|T§|l becomes: 



q n\ __ e i6x ~ A - e i6 £ 

k\ TIq e M\-A _|_ e iS t 



where we have defined 



= tan 



Ox-A - i 



, (20) 



9\~a = arg(A — A) , 6 S = arg(B). (21) 
To determine the latter, we note that 



n\k 2 



n 2 ki 
n\k 2 



(22) 



appearing in eq. Q can have either sign. In fact we 
expect that the sign will change when /3 runs from tiq = 1 
to n 2 . We therefore write 



24) + S 



(23) 



with </> s chosen to be or — 7r depending on the sign of 
c_ sin(k 2 d 2 ). As shown in Ref. 10], this sign is char- 
acteristic of each bandgap, unless the width shrinks to 
zero in either an optical hole [Tol ] or a Brewster point. 
Returning to eq. I|20|l . we can now write 9\-a explicitly 
as 

X -A(k) = 20(/3) + kmpfa - 24) + cb s + | + 2tt^ (24) 
with v = 0, ±1, ±2, ... , and where we have defined 



III. SURFACE TM WAVES 

Most derivations are analogous to TE case, see Ref. 
:9j. The magnetic field of a surface TM wave is written 
as 

H y (z) = a s e^ z +b s e- tk ' z , (17) 
when —4 < z < and 

H y {z) = b v e qoZ , (18) 



with qo = +ky/ j3 2 — tiq, when z < —4- Using the 
boundary conditions, we obtain the exact dispersion re- 
lation k = k((3) for TM surface waves by solving 



qo_ n? 
B™ 



i = . 

-2ikid s ^jTM 



(19) 



This equation must be solved numerically, and we will 
refer to the solutions thereby obtained as "exact". We 



Q(f3) = taa 



ll3 2 -nl 



(25) 



Note that all dependence on k on the r.h.s. of eq. Ij24(l 
is explicit in the second term. Note also that in taking 
the inverse tangent we have introduced a contribution of 
+27rf on the r.h.s. This term was discussed at length in 
Ref. 0, so here we will discuss only the solutions with 
v = 0. 

There is no simple analytic form for Q\~a as a function 
of A;, so numerical methods must be used to solve eq. I|24ll . 
However, the results can be better understood using a 
simple graphical approach. In Fig. [21 we plot 9\-A(k) 
(continuous line), and the r.h.s. of eq. I|24|l for a chosen 
set of 4's (dashed lines). One sees that 9\-a increases 
from — 7r/2 to n/2 as k ranges from the lower to the upper 
first bandgap boundaries. In a more general context, this 
property has been shown to hold in Ref. ^1 • For a given 
4, the intersection of the corresponding straight line with 
the continuous line determines the solution for k. As the 
graph shows, when 4 decreases, k increases. The values 
of 4 for which a solution can be found will therefore be 



bounded by those corresponding to the r.h.s. of eq. i|24|) 

where it crosses Q\-a — ±7r/2. That condition leads to 



beta = 1 .01 



1 



di 

2 k R nij 3 

d x 1 
— H 

2 k L ni t p 



e(/3) + y ) i 



e(/3) 



2 2 



(26) 



In the example shown in Fig. |2 one finds d c ,min = 21 
nm and d c ^ max = 134 nm (For this example, cj> s = —tt 
and (3 > 0Br)- 



beta = 1 .30 




k ( nm -1 




6^067 O.O'OSS 0.0'069 0.607 0.0'071 0.0'072 0.0073 
k (nm-1 ) 



FIG. 3: Graphical solution of eq. I|24|l for no = 1, n\ = 
2, di = 100 nm, n 2 = 1.5, d 2 = 250 nm, and (3 = 1.01. 
Continuous line: 8\-a', dashed lines: r.h.s. for several values, 
from top to bottom, d c = 0.90, 1.20, 1.70 and 2.0di. 



uq to rt2- For this example, the approximation is so good 
that on that scale one cannot distinguish between the 
exact and the approximate curves. 



FIG. 2: Graphical solution of eq. I24H for no = 1, ni = 
2, di = 100 nm, n 2 = 1.5, d 2 = 250 nm, and j3 = 1.3. 
Continuous line: 8\-a; dashed lines: r.h.s. for several values, 
from top to bottom, d c = 0.25, 0.90, 0.99, 1.01 and 1.33di . 

The same expressions also apply when p < /3s r , but 
then (b s = 0. If we choose (3 = 1.01, the limits are 



75 nm and d c 



134 nm. The graphical 



solutions are now shown in Fig. [3] 

To proceed further in the analysis of the SEW solutions 
requires values for the band edges fci, kp> for a given /3. 
These can be obtained from our semiclassical approxi- 
mation [§[■ In Ref. _9j we found very accurate analytic 
approximations for the TE bandgaps. Formally analo- 
gous expressions hold for TM bandgap boundaries, and 
their accuracy is also excellent. We have also shown that 
the semiclassical theory allows one to derive a good ap- 
proximate expression for the argument of A — A. In the 
first bandgap, it is 



a(e) _ 

1 \-A — 



sin 



w/2 



(27) 



with k m = (kpi + kL)/2 and w = kn — k^. 

In Fig. 2] we compare the exact valand the above ap- 
proximation for Q\~a for several values of f3 ranging from 




FIG. 4: Exact (continuous line) v.s. empirical approximation 
(dashed line) to the argument of A — A. From left to right, 
f3 = 1.05,1.1,1.15,1.25,1.3,1.35 and 1.4. We omitted the 
curves corresponding to the Brewster point at 0Br = 1-2 
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A. Approximate analytical solutions in the middle 
of the first bandgap. 

SEW in the middle of a bandgap are of interest because 
the damping is strong, so most of the surface wave is 
confined very close to the surface. For momenta such that 
k — k m << w one can make the simplification sin -1 (2 (A; — 
k m )/w) ~ 2[k — k m )/w and eq. j^H then gives 



6(/3) + k m /w + <f> s /2 + ttQ/ + 1/4) 
1/w + n lj3 (d c - tii/2) 



(28) 



The dependence of k on d c , (3 and d\ is now transparent. 
The role of the bandgap parameters k m and w can also 
be easily discussed. Fig. [3] shows the accuracy of this 
approximation when (3 > (3bv 




0.022 
0.02 
0.018 
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FIG. 6: Dispersion relation k — k((3) when d c — 0.60 di . Ex- 
act solution (continuous line) v.s. approximations : Crosses: 
eq. il'i.'-il : stars: eq. $i4) . The dashed line is the upper bound- 
ary of the first bandgap. 



Inserting this into eq. (|24|l gives 

2 jkR-k + A 
V w 

with 

p = n ip {2d c - di) A = -26 - cp s + k R p. (32) 
Solving for kjj — k, we obtain 



FIG. 5: Dispersion relation k = fc(/3) when d c = 0.75 di, 
0.90 di. Continuous line: exact solution of eq. 1191 . Dotted 
line: linear approximation, eq. 1281 . Crosses: first bandgap 
boundaries. 



B. Approximate solutions near the bandgap 
boundaries. 



Here we discuss solutions near the upper bandgap 
boundary, but similar approximations can be developed 
for the lower boundary, as in Ref . Q . When k is slightly 
below kg. = k m + w/2, it is convenient to introduce £ > 
via 



fe m = -(i-C). 



(29) 



k = k R -^l-^=+\l-+ P .\ 



if 



(33) 



which is the desired solution, k = fc(/3), near the upper 
bandgap boundary. Fig. shows an example of the ac- 
curacy of this expression. Furthermore, when pA is small 
compared to 1 /w one can expand and find 



k = kj 



-wA 2 



(34) 



which again manifests the dependence of k on d c and w, 
and allows one to construct k — k(f3) very easily. Fig. 
shows again the validity of this approximation. 



Then, 



sm 



. w/2 



- 
2 



■2, = l-2J k -^. (30) 
2 V w 
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IV. FIRST AND SECOND SEMICLASSICAL 
APPROXIMATIONS 

The transfer matrix relates the amplitudes for trans- 
mission t B and reflection r B of a single cell as 



M 



1 

TB 



(35) 



Comparing to eq. {J, A = l/t* B and B = —r B /t* B . In 
the first approximation, the expressions for r B and t B 
are given in Refs. 0,0] for TE waves and in Ref. 8] for 
TM waves. In particular, we have 



-s* sinh(7id) 



7i cosh(7id) — iS q sinh(7itf) 

tYhi 

71 cosh(7id) — iSq sinh(7id) 



(36) 



where the coefficients s q differ of course for TE and TM 
waves. For a photonic crystal with a symmetric cell, 
d\ = 2dt, they take the form 



„TE 



In 



sin 



[qn + (k 2 - ki)di 



(37) 



„TM 



= --.In 



111 


ki n 2 


e- tq7r sin 




k 2 n\_ 





-J [?7T + (fc 2 



h)di 



The parameter 5 q is the detuning from the q-ih Bragg 
resonance k q , q = 1,2... and the well-known Bragg con- 
dition for constructive interference is 



q, S q = kn av 



(39) 



The parameter 71 appearing in eqs. (|36() is defined as 
71 = \J\s q \ 2 — S 2 . Then, for the elements of the trans- 
fer matrix in the first approximation of the scmiclassical 
coupled wave theory, we obtain 



A {1) = (-)« (cosh7id + i^sinh-,(/ 



B (i) = (_) 9 fi sinh7ld , 
7i 



(40) 
(41) 



Inserting Re (A^) into eq. (fTT|) . we find that A^ = 
(— ) q e~ lld , where A' 1 ) is the first semiclassical approxi- 
mation to the eigenvalue of the damped Bloch wave. 



The second approximation introduced in Refs. Q 
leads to similar expressions for the elements of the trans- 
fer matrix: 



A< 2 



cosh 72 d 



+ 1 



[(1 + M 2 )^ - 2lm( Sq u*)] 
(1-M 2 )72 



sinh72<i 



B m 

where 72 = 



s q - 2ir\ q u - s*u 2 . 

r-r^ smh 72 d , 

(1 - \u 2 72 



,(42) 



(43) 



s <?l 2 _ Vqi Vq = Sq — ic 2 and 
id 



C2 



E 



27r ^— ' m — q — Snd/Tr' 

■m^q— — 00 



m— +C3Q 



7T 



2n ^ m — q — Sgd/ir 

m= — oo, m^q 



(44) 



(45) 



The eigenvalue of the damped wave is A^ 2 ^ = (— ) q e l2d . 

Inserting A^ 1 ' 2 ) and B^ 1,2 ^ into eqs. (|19|l . one finds the 
corresponding predictions for k = k(f3) of TM surface 
waves. Selected results are shown in Fig. As can be 
seen, the second approximation gives excellent agreement 
with the exact dispersion relation within the range of 
validity (3 < n 2 of the semiclassical coupled wave theory, 
for the problem at hand. 
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FIG. 7: k — k(/3) curves for TM waves in the first bandgap 
for d c — 50, 75 and 90 nm; continuous lines: exact solutions; 
long dashes: first semiclassical approximations; short dashes: 
second semiclassical approximation; crosses: bandgap bound- 
aries 



7 



0.022 




As shown in Ref. 0, the dispersion relation k = k((3) 
for TE surface waves is found by solving 

^ ^ A TE „ A TE „ TE 

fa ~ ~ l A™ - A TE + B TE ' 
B TE = e - 2ifeld 'S TE . (46) 

In Fig. El we compare the exact and semiclassical results 
for k = k(/3), choosing three thicknesses of the cap layer 
d c = 25, 50 and 75 nm. The first approximation becomes 
inaccurate when (3 exceeds 1.4, but gives accurate results 
up to that value. The second approximation is so close to 
the exact values that one can see the difference only for 
values of /3 very close to the critical value j3 = 712 = 1.5. 
Beyond that, our semiclassical approximation cannot be 
applied, since /C2 becomes imaginary. 



FIG. 8: k — k(/3) curves for TE waves in the first bandgap 
for d c = 25, 50 and 75 nm; continuous lines: exact solutions; 
long dashes: first semiclassical approximations; short dashes: 
second semiclassical approximation. The latter curves are so 
close to the exact ones that the difference can be seen only 
when (3 > 1.45. 
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